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Abstract 

In the context of spin electronics, the two spin-channel model assumes that the spin carriers are 
composed of two distinct populations: the conduction electrons of spin up, and the conduction 
electrons of spin down. In order to distinguish the paramagnetic and ferromagnetic contributions 
in spin injection, we describe the current injection with four channels : the two spin populations 
of the conduction bands {s or paramagnetic) and the two spin populations of the more correlated 
electrons {d or ferromagnetic). The redistribution of the conduction electrons at the interface is 
described by relaxation mechanisms between the channels. Providing that the d majority-spin 
band is frozen, s — d relaxation essentially concerns the minority-spin channels. Accordingly, 
even in the abscence of spin-flip scattering (i.e. without standard spin-accumulation or giant 
magnetoresistance) , the s — (i relaxation leads to a d spin accumulation effect. The coupled diffusion 
equations for the two relaxation processes {s — d and spin-flip) are derived. The link with the 
ferromagnetic order parameter M is performed by assuming that only the d channel contributes 
to the Landau-Lifshitz-Gilbert equation. The effect of magnetization reversal induced by spin 
injection is explained by these relaxations under the assumption that the spins of the conduction 
electrons act as environmental degrees of freedom on the magnetization. 

PACS numbers: PACS numbers: 75.40. Gb, 75.70. Pa, 75. 70. Cu 



* Electronic address: |jean-eric.wegrowe@polytechnique.fr| 



1 



Spin-dependent transport in intermetallic ferromagnetic materials is usually described 
in the framework of the two spin-channel approximation. In this model the spin carriers 
are composed of two distinct populations: the conduction electrons of spin up (t), and the 
conduction electrons of spin down (|), which are defined by their conduction properties 
This model can be applied to ferromagnetic junctions and to hetero-j unctions and 
it provides simple and operational descriptions of many effects related to spin dependent 



transport. The giant magnetoresistance (GMR) effect j2| 



IB 



is described by the 



diffusive spin-accumulation mechanism in normal / ferromagnetic junctions 



and heterojunctions 



I2I ] and also by interface scattering and/or band mismatch effects. 



Beyond the well-known GMR, which is due to the action of the magnetization on the 
spin of the conduction electrons, the opposite effect of the spin polarized current on the 

ed "current induced magnetization switching" or CIMS jisl 



ma gne tization, the so caf 
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22] has also been treated in the framework of the two channel 



approximation (see the refrences [2^ for the effects related to noncoUinear magnetization). 
Tremendous attention has been paid to the microscopic description of CIMS in terms 



of the dynamics of the spin of the conduction electrons in ferromagnets 



29 



^12^12 



\21 



|. These approaches add a deterministic correction ("the spin torque") to the 
Landau-Lifshitz-Gilbert {LLC) equation due to the current, by identifying the spins of 
conduction electrons to the ferromagnetic order parameter. The old and unconfortable 
discussion about paramagnetic vs. ferromagnetic character of itinerant spins (already 
problematic at equilibrium as discused in the seventeen j^ll) is implicitly evacuated. A 



care: 



uU study, however, of the different kinds of experiments performed about CIMS 
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38, 



39 



4G 



41 



42 



|45| motivated in parallel the devel- 



oppment of a stochastic approach, where the action of the spins of the conduction electrons 
on the magnetization is not direct, and is described in terms of environmental degrees of 
freedom. The magnetic order parameter and the spins of conduction electrons have not the 
same physical meaning. The most typical stochastic effect is the two level fluctuation and 
activation processes observed only under current injection over a time window range of a 



few nanoseconds 



35 



38, 



to minutes jl^. Below the nanosecond time range, time 
resolved measurements (i.e. one-shot measurements) are not available but absorbtion peaks 
are measured in the electric power-spectrum at frequencies between 3 to 20 GHz 0, |3| ■ 
These resonances can be though of as a reminiscence of what is observed in ferromagnetic 
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resonance experiments, or noise measurements without spin-injection :at these time 

scales only (sub-nanosecond), thermal spin- waves and precession of the magnetization are 
expected. This is the typical time scale of the dynamics of the magnetization. At the smaller 
coarse-grained level, electronic relaxations take place at times scales ranging from 10^^^ to 
10~^^ second. These well separate time scales justify an approach in terms of Markovian 
processes in which the relevant slow variable is the ferromagnetic order parameter, and the 
fluctuations come from the relaxation of electronic spin degree of freedom. 

The aim of this paper is to propose a phenomenological description of a possible mech- 
anism of electronic relaxation (presented in Section I) which accounts for the fondamental 
difference between paramagnetic and ferromagnetic spin currents. In contrast to the para- 
magnetic current associated to the standard spin accumulation process and GMR, which 
does not interact direcly with the magnetization, the ferromagnetic, or d channel current 
(the denomination "d channel" is a reference to the equilibrium description of itinerant d 
electrons 0|) contributes to the ferromagnetic order parameter. The relaxation from 
paramagnetic channels {s channels) to the ferromagnetic channels {d channels), described 
as a diffusion process at the interfaces, leads to spin-transfer from injected current to the 
d current, through a new d spin accumulation process (Sections 1 and 11). It is furthermore 
assumed that even this d contribution to the dynamics of the magnetization is not direct, 
due to the difference in the typical time scales. The d current injection is accounted for in 
terms of magnetic fluctuations or noise (Section 111) and the after effect is to excite a large 
spectrum of magnetic and non-magnetic modes in the ferromagnetic layer. 

This approach accounts for both the measured stochastic processes of the magnetization 
with a cut-off frequency beyond one GHz, and the observation of magnetization reversal by 
spin injection measured in sing le uniform magnetic layers (i.e. in the abscence of giant mag- 



netoresistance) 



38, 



39I |40| . The approach adopted here is based on general arguments 



: 'rom irreversible thermodynamics, following the general method of De Groot and Mazure 
3 or Stuekelberg 0. 

Following the pioneering ideas of Mott [53], we assume that two electronic populations, 
referred to as s and d are relevant in order to describe the conduction properties of the 
ferromagnetic 3d metals and alloys. This hypothesis was used in order to describe the 
anisotropic magnetoresistance in these compounds 1541, and to account for the contribution 

nn 

to the ferromagnetism of itinerant d electrons in equilibrium |42l, |^ . In our context, the 



3 



s population describes the conduction electrons injected at the interface of a ferromagnetic 
nanostructure, and the d population describes the more localized electrons and associated 
spins, mainly responsible for the ferromagnetism of d ferromagnets. The simplest extension 
of the two-channel approximation to our context is to divide both the s and d populations 
into up (t) and down (|) spin populations (| is the internal variable describing the spin 
degree of freedom as defined in Ref. 5^). 

Providing that the d majority spin band (|) is fully occupied, s — d relaxation of | spin 
is negligible, and spin-conserved s — d relaxation only concerns the minority spin channels 
|. If we assume further that s — d relaxation with spin-flip is negligible with respect to spin 
conserved s — d relaxation, a new mechanism of spin injection could then be described as 
parallel to the usual t — i spin-flip relaxation used to describe spin accumulation and GMR 
phenomena 0, 0, [l^ . 



I. THERMOKINETIC EQUATIONS 

The system is composed by the reservoirs of the injected s electrons and the ferromag- 
netic layer composed by the d electrons. At the interface, current injection leads to a 
redistribution of the different electronic populations that are governed by spin polarization 
and charge conservation laws. Let us assume that the current injected is spin polarized in 
the down polarization (|). The conservation laws should be written by taking into account 
the reaction mechanisms between the different populations. At short time scales (electronic 
scattering) the relaxation channels are assumed to be the following four 



(I) — > Cdi (spin-conserved s-d scattering) 

(II) (spin-flip scattering for the s population) 

(III) 6(^1 (spin-flip s-d scattering) 

(IV) (spin-flip scattering for the d population) 



In a second step (at larger time scales), opposite relaxation mechanisms take place (in- 
volving ferromagnetic excitations) that are described in Sec. III. Process (I) is assumed to 
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be the main mechanism responsible for irreversible spin transfer. Process (II) leads to the 
well-know n sp in-accumulation effect and was described in detail with the same formalism 



elsewhere jl5| . According to the fact that the majority- spin d hand is full and lies at a sizable 
energy below the Fermi level, the current Jd] i s negligible and the channel (i f is frozen. 
Processes (III) and (IV) are hence negligible [56[ . Consequently, we are dealing with a three 
channel approximation. 

The total current Jt is composed by the three currents for each channel {s t, s [,d [}: 
Jt = Jsl + Jsl + Jdl- In order to write the conservation laws, the relaxation rate "^sd, is 
introduced to account for s — d spin-conserved scattering, and the relaxation rate is 
introduced in order to account for spin-flip scattering. Assuming that the current is flowing 
along the z axis, the conservation laws for a steady state regime are : 



dJt _ 


dJt _ 





dt 


dz 


dJal 

dt ~ 


dJs^ 
dz 


- = 


dJsl 

dt 


dJsi 
dz 


- + = 


' _ 

dt 


dJdi 
dz 


+ ^sd = 



(1) 



The system is described by the number of electrons present in each channel at a given time, 
that defines the four currents, plus the entropy of the system. The conjugated (intensive) 
variables are the chemical potentials {^si, fJ^si, fJ'di, fJ^di}- As described in Appendix A, the 
application of the first and second laws of thermodynamics allows us to deduce the Onsager 
relations of the system : 



J 



*i e dz 

"'^T — e dz 

U = (2) 

^sd = Lsd ifJ'sl — fJ'dl) 
= Ls if^sl - f^sl) 

where the conductivity of each channel {asi,<Jsi,o'di,o'di} has been introduced. The 
first four equations are nothing but Ohm's law applied to each channel, and the two last 
equations introduce new Onsager transport coefficients (see Appendix A), Lgdi and Lg, that 
respectively describe the s — d relaxation (I) for minority spins under the action of the 
chemical potential difference A/xj^ = ^si — fJ-di and the spin-flip relaxation (II) under spin 
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pumping A/is = M - l^si- 

The quantities of physical interest are the paramagnetic current Jqs = Js^ + Jsi, the 
minority-spin current Jqi = Jsi + Jdi, and the two polarized currents 53^ = Jsi — Jdi and 
5Js = Js-\ — Jsi- We introduce the ag and o"| conductivities {ag = ag-^+cXsi and = o'si+cTdi}- 
The conductivity imbalance Pi and Ps between respectively the s | and d | channels and 
the s t and s | channels are: 



(3) 



Eqs. (HJ becomes 



dJt 



dz dz 8. 



+ 



-2ijsd - (4) 



asjy 

dz ~' 
'-t = -^sd 

^ = - 2^, 

and, defining the quasi-chemical potentials |3| A^s = /^st + /^«i ^^"^ A^i ~ /^si + A^rfi; Eqs. (j2j) 
becomes : 

= -ft {Ps'-t + ^) 

^fsd = LsdAfil 

The equations of conservation [Eqs. and the above Onsager equations lead to the two 
coupled diffusion equations : 

SO. sj 
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where 



- Y 2eL, 



(7) 



4eL, 
gs(l-fe) 



A solution of Eqs. © is 



(8) 



with 



A/ii = aie'^+ +026 
A/i2 = he^ + &2e 



A_ 



(9) 



where 



Ai^ = ^(C.^ + C/)|l± 



/-2/-2 _ \-2\-2 



The constants ai, 02, &i, &2 are defined by the boundary conditions. It can then be seen 
that the usual spin accumulation corresponding to A/i^ also depends on the spin-conserved 
s — d electronic diffusion which is known to be efficient and, conversely, that spin- 
conserved diffusion is able to lead to a spin accumulation, or d spin- accumulation effects. 
Accordingly, we expect to measure some typical effects related to spin-accumulation in single 
magnetic layers, or if /^^ = : this point will be illustrated in the new expression of the 
magnetoresistance (Eq. (fT^ below), and in Section III through the effect of CIMS. s — d 
relaxation adds a new contribution to the resistance, which plays the role of an interface 
resistance arising from the diffusive treatment of the band mismatch 

The resistance produced by the spin-ffip contribution (usually defined as the giant mag- 
netoresistance Rgmr), plus the contribution of s — d relaxation, are defined by 



7 



where Et = — ^ is the total electric field and E'j^ is the electric field far away from the 
interfaces. Providing that the total current is Jt = Js^ + Jsi + Jdi, or 

(^t d f (Tsi 0-4 



e dz \ at (Jt at ' J 

The total electric field can also be written (from Eqs. Q) as 

= = :Ii = -- (^t^ + ^^^^^ (12) 

edz at e \ dz at dz at dz J 

where at = agf + agi + aai , and = limz^^oo ~ The resistance is given by : 

B_._J_p(-^ + !-l^?§^^i),, (13) 

eJt J -00 \crt dz at dz J 
This three-channel model brings to light the interplay between band mismatch effects and 

spin accumulation, in a diffusive approach. It is interesting to note that the local neutrality 

charge condition which is often used (see for instance Eq. (4) in |5S]) was not included. On 

the contrary, we have imposed the conservation of the current at any point of the conductor. 

Indeed, electron transfer from a channel to another where the electron mobility is different, 

induces a local variation of the total current. 



II. NORMAL - FERROMAGNETIC INTERFACE 

The resolution of Eqs. (jHI) leads to a variety of possible behaviours, from single inter- 
face effects to superlattice effects (see the paper by Valet -Pert for the discussion in the 
framework of the two-channel approximation). Our main goal however is to understand the 
contributions of standard spin-accumulation and d spin-accumulation in uniform magnetic 
layers where no GMR are present due to the symmetry between both interfaces. We first 
focus our attention on a single interface separating two semi-infinite layers. If the ferromag- 
netic layer is at the right hand side, the solutions follow, from Eqs. (jH)) 

A^i(z) = be'^ ^^^^ 
Afi2{z) = de~^ 

In our context, we define a "normal" metal as a compound with fully occupied d bands 
(the d I and d | channels are frozen). In the normal metal (on the left hand side of the 
junction) we have Afi^ = ae~^/'=^. 
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At the junction, the continuity of the currents and the continuity of A/i^, without interface 
resistance, becomes : 







^il(o-) = ^fi(0+) (15) 
A/if (0-) = A/if (0+) 

Where the superscripts N and F resp. stand for normal and ferromagnetic metals. The 
system originated from Eqs. ()15|) is solved in Appendix B. The three limiting cases of 
special interest are presented below. 

Let use first consider the standard assumption where GMR and spin-accumulation are 
calculated. This case corresponds to both limits Isd <^ hf (i- e. instantaneous s—d relaxation 
at the interface), or l^f <^ /^^^ (no sd relaxation in the time scale of spin-flip relaxation). The 
simple case of two identical adjacent layers with opposed spin polarizations is generally 
considered P|. In our context, this situation would correspond to a " quasi- ferromagnetic" 
metal in which both d | and d | are frozen (i. e. there is no d ferromagnetism) but with 
/9f = —(3s^i i. e. there is nevertheless a spin polarization of the current. The well-known 
solution is straightforwardly recovered. From the equations presented in Appendix B: 



2e(3,ls_ 

In the case of a junction with a "normal" metal /3f = and a ferromagnetic metal jSg 7^ 0, 
the well-known result is also recovered in both limits (with erf = = at and Ifj- ~ Ifj-) : 

A/^s - ^^^^*e-^/W (17) 
(1 - ^) 

in these two limits there is hence no modification due to the existence of the "d spin- 
accumulation" of the standard spin accumulation A/i^. The d spin-accumulation however 
is not zero. In the case of a normal metal / ferromagnetic junction with , the d 

spin-accumulation is : 



The corresponding contribution to the magnetoresistance is proportionnal to l^^ and 



— W 

^ 2 



This contribution is not zero even if the spin polarization of the current is vanishing (/3s = 0) 



in the ferromagnet. Since Igd is however expected to be very small, the d spin- accumulation 
contribution should be important only for Ps close to (this is a consequence of the 
assumption cr^lff = a'^lff). 

In intermediate cases, when Isf is of the same order of magnitude as Isd, the spin- 
accumulation Afis is non-vanishing even if (3s tends to zero. 



III. GENERALIZED LANDAU-LIFSHITZ-GILBERT EQUATION 

The question related to magnetization reversal under spin-injection is how the kinetics 
of the spins of the conduction electrons described in the previous sections is related to the 
ferromagnetic order parameter. In order to investigate this problem, we shall consider the 
different quantities that define our system (spins of conduction electron, electric charges, 
magnetization . . . ), and the corresponding relevant scales, or coarse-graining 

In our experimental context, the magnetization is a collective variable whose dynamics 
are much slower than all other paramagnetic spin relaxation mechanisms because the mag- 
netization is conserved over the distance of the magnetic layer (or exchange length for spin 
waves), while the conservation of the electric charges is relevant over a local equilibrium 
of the order of the nanometer. It is then possible to identify three different time-scales. 
The first is the electronic relaxation (10~^^ to 10~^^ for paramagnetic transvers spin effects 
[2^), the second is the typical ("quasi-ballistic") dynamics of the magnetization (10~^^ to 
10~^), and finally, the time scale of the activation processes measured over decades from 10 
nanoseconds to hours. The electronic degrees of freedom, and especially the spins of the four 
electron populations defined in the first section can then be treated in terms of the action 
of an environment in a stochastic approach (e. g. defining a projection operator over the 
relevant variable Mq The effects of the spins of the conduction electrons are 

then reduced to the noise and the damping coefficient. Without current injection Q], the 
electronic degrees of freedom are contained in the Gilbert damping term |62l . 163^ , and the 



bl 



owing Gilbert equation for the magnetization Mq of the ferromagnetic layer is obtained 



dMo ^ dV dMo[ 

— — = r Mo X < ri —— > (19) 

dt \ dMo dt j ^ ^ 

where T is the gyromagnetic ratio and V is the magnetic Gibbs potential j6^], where 
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^ext ijging |;]2e external magnetic field. The r] coefficient is the Gilbert damping factor. 
Eq. ()19|) can also be put into the following Landau-Lifshitz form. In the case of uniform 
magnetization we have Mq = MgUQ, where Ms is the saturation magnetization and Uq the 
unit vector, so that it becomes: 

^ = - ^' (^uo X VV^ - h'uo X (uo X VV^ (20) 

where V is the gradient operator on the surface of a unit sphere. The phenomenologi- 
cal parameters h' and g' are linked to the gyromagnetic ratio F and the Gilbert damping 
coefficient r] by the relations 



h' 

9' 

a 



Fa 

r 

(l+a2)M, 



If we take into account the relaxation processes (I) to (IV), the conservation of the 
ferromagnetic spins normal/ferromagnetic interface is: 



u 



g' (ux VV) - h'u X (ux Vv) - fii{t)u 



(21) 



where 



9dP'B 



((5Jf *(2) - 5Jf"*(2)) dz 



(22) 



L is the length of the magnetic layer, and QdfJ'B is the magnetization per atom (number of 
Bohr magnetons) for the d population; gfis is the magnetization per atom of the ferromagnet, 
and the superscripts int and out describe respectively the interface for the incoming current, 
and the interface corresponding to the current flowing out of the layer. The (5Jj current is 
defined in Eq.((Hl) : 



6Ji 



2e dz dz 



(23) 



with the solutions given in Section II and Appendix B. Note that u is no longer a unit 
vector due to the last term in the right hand side of the equation Eq. Instead, u includes 
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the fluctuations of uq. After projection over u, 
of the magnetization due to current injection ( 



Eg . (PT|) gives the variation of the modulus 
3): 



d\\uf{t) 



dt 



-— llniP 



(5Jf - SJl'^'iz)) dz 



(24) 



and 



\\u\\\t) = ||«oire-{fi^o"('^r*(^)-'^r*W)'^4* (25) 

The integral in the exponential is rather similar to that present in the calculation 
of the GMR, except that the current is the ferromagnic spin current. For symmetry 
reasons equal and opposite electronic relaxations are expected to occur at both inter- 
faces of the ferromagnetic layer : the effect of spin-injection should be compensated : 
/ i^Jl^^i^) ~ / SJ°'^^{z)^ dz = 0. We may still conclude that, if d [ spins do not interact 
with the magnetization and if the two interfaces of the ferromagnet are totally symmetric, 
the spin-transfer would vanish. 

This is indeed what happends in the case of s population and GMR, where the energy 
due to the spin-flip is tranferred into the lattice and dissipated into the heat bath. This is 
simply due to the thermalization associated to paramagnetism. To that respect, the idea of 
conservation of momenta which leads to spin-transfer theories is not valide for paramagnetic 
spins, because the momentum dissipates into the heat bath and is not transfered into the 
ferromagnetic order parameter. In Eq. 1)2111 . the transvers stochastic force responsible for 
the fluctuations f±{t) has been put to zero because of the averaging over the equilibrium 
distribution {{f±{t))o = 0): the transvers stochastic force is suppose to be the same with or 
without current injection and (/±(0)/_L(it))o = crjkTdit) where T is the temperature of the 
lattice and c = VMg 

In contrast, the energy dissipated or gained hy s — d scattering is tranfered in the form 
of a contribution to the ferromagnetic order parameter within a typical time scales of the 
dynamics of the magnetization, and is not dissipated into the heat bath. Thus, the d [ 
spins do interact with the magnetization within a short distance. The process described 
in Eq. ()^3|) can only be observed locally, and in the sub-nanosecond time range. Beyond, 
the behavour of the magnetic layer will reflect the large spectrum of relaxation channels 
of the magnetization, from spin-waves, solitons, precession, soft modes etc. The detailed 
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description of these contributions is beyond the scope of the present paper, but the resulting 
behaviour of the magnetization can be described phenomenologically in terns of random 
fluctuations of the magnetization and effective temperature. Instead of Eq. (j25|) . we have 
= limt^oo II^P(^) ~ 0. At long enough time scales or thick enough layers, the last 
term in the right hand side of Eq. (j2ip is averaged out : 

(/iiW)^o 

The stationary transfer of spins, described here as environmental degrees of freedom, 
emerges in terms of transfer of energy (and entropy) through the non-vanishing fluctuations 
of the modulus of the magnetization. The fluctuation- dissipation relation, defines the current 
dependent effective temperature Teff{Jt) with the correlations: 



{f{0)\\f{t)^\)^cvkT,ff{Jt)S{t) (26) 

where k is the Boltzmann constant and c is an appropriate constant. The energy kT^ff 
is stored in the layer in form of magnetic fluctuations, and is not direcly related to the 
temperature of the lattice. Accordingly, it is possible to transfer some few eV without 
damage in a nanoscop.c sample (except .f the fiuctuatiou^are generating spin-waves only, 
because spin- waves relax very rapidly into the lattice) [32]. It is also expected that the 
efficiency of the transfer is maximum if the layer thickness is large. This would explain 



the 



IS 



Dehaviour as the function of the temperature observed in ultrathin trilayer nanopillars 
3q : the amplitude of CI MS decreases while decreasing the temperature, and in long 



Ni nanowires [39| where the amplitude of CIMS is constant or increases while decreasing 
the temperature. 

The effective temperature is proportional to , and depends, through A/i^ on the spin 
accumulation properties, and through AyU^ to the d spin-accumulation. The fluctuations 
of the current SJi{t) are hence able to account for the effective temperature measured in 



some experiments of magnetization reversal 32, |3J, |35|, |38[ in both single layer and trilayer 



systems. The effect described in terms of effective temperature is not directly sensitive to 
the sign of the current, but depends indirectly on the current direction through the sign of 
A/is and A/i|, because it inverses the asymetry at the interfaces jsj in multilayered GMR 
devices. In contrast, the spin accumulation is symetric in a single magnetic layer so that the 
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asymmetry of the interfaces is not necessarily modified by inverting the current. It is then 
)ossible to account for both the dependence to the current direction in trilayer structures 

19, 2^0,0 QQ, 

and the absence of dependence in AMR single layer measurements 

37, 




On the other hand, the resonance measured at the GHz frequency range j36l . l4lL |42| | may 
be described by the magnetic excitations (precession and spin-waves, etc) produced during 
the relaxation of the modulus of the longitudinal magnetization u into the ferromagnetic 
layer, while direct measurements of electronic resonances should be expected at the 10 to 
100 GHz frequency range. 

It is not possible to perform a direct estimation of spin transfer at the present stage of the 
description because the amplitude of the fluctuations strongly depends on specific magnetic 
properties of the sample and interfaces. A rough estimate can nevertheless be performed by 
observing that the typical current density injected at about 1.6 lO'^A/cm'^, which corresponds 
to 10^^ electrons per second injected at the interface. The current Joi is a fraction of that 
(lets say above 1 % The question is to know what is the typical length 1, or the 

time T (l=v r ) over which the magnetization is maintained out of equilibrium in the d | 
channel. The minimum relaxation time should be about 10^^^ seconds, the maximum should 
be the typical magnetization dynamics, around 10"^ seconds. This means that the system 
is pertubated by a magnetization variation of 10^ to lO^/i^. Note that 2.10'^ fis corresponds 
to an energy transferred by the current of more than one eV in a local ferromagnetic field 
of 1 Tesla (in agreement with experimental results js^). It is then possible to account for 
a transfer of magnetic momentum, with an energy largely beyond the Curie energy in the 
local field, and which would produce the magnetization reversal or magnetic excitations. 



IV. CONCLUSION 



A new electronic relaxation mechanism has been proposed under spin-injection at a nor- 
mal/ferromagnetic interface. The description of the relaxation is based on a three-channel 
model that leads to a redistributions of electrons between paramagnetic and ferromagnetic 
spin currents at the interfaces. As a consequence, a new spin accumulation process of the d 
electrons occurs. The coupled diffusion equations are derived and solved. The contribution 
of the d spin accumulation to the standard spin accumulation in the GMR is calculated. 
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The d spin accumulation adds a new contribution to the resistance, which plays the role of 
the interface resistance arrising from the diffusive treatment of the band mismatch. 

In contrast to the paramagnetic current (here associated to the standard spin accumula- 
tion process) which does not interact direcly with the magnetization, the ferromagnetic, or 
(i-channel current contributes to the ferromagnetic order parameter. It is furthermore as- 
sumed that even this d contribution to the ferromagnetic order parameter is not direct, due 
to the difference in the typical time scales. The d current injection is accounted for in terms 
of magnetic fluctuations or noise the consequence of which is to excite a large spectrum of 
magnetic and non-magnetic excitations in the ferromagnetic layer. This stochastic approach 
allows an effective temperature (or equivalently an effective potential barrier) to be defined 
in agreement with the experimental observations. The fluctuations depends, through the 
interband current 6Ji, to both the usual spin accumulation Afig = fisi — /^sT ^ 
spin accumulation A/i| = /i^i — /irf^. This mechanism allows the effect of current induced 
magnetization switching, including current-induce activation, to be described not only in 
multilayered structures exhibiting GMR, but also in uniformly magnetized nanostructures 
measured with AMR. 
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VI. APPENDIX A 

The aim of this Appendix is to derive the Onsager matrix (0) on the basis of the first 
and second laws of thermodynamics. In a typical one dimensional junction the layer is 
decomposed into Q parts, defining the sub-system "E^, which is in contact to the "reservoirs" 
and The sub-systems S'^, is then an open system which exchanges heat and 

chemical species with its left and right vicinity layers. Furthermore, the populations (A^^|, 
N^^, N^.^) and spin down (N^^) are not conserved due to transitions from one channel to the 
other. 
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In this picture, the states of the sub-system E'^ are described by the variables 



{S\N^^,N^^,Nl^,Nl^) (27) 

where 5"*^ is the entropy. The internal variables and must however be intro- 

duced in order to take into account the relaxation processes occuring respectively between 
the two s-spin channels, the two d-spin channels, and the s-d relaxation . 

Let us define the heat and chemical power by (the mechanical power is zero as long 
as the action of the magnetic field on the charge carriers is neglected). The first law of the 
thermodynamics applied to the layer E*^ gives 

Introducing the canonical definitions = and = -§§w--> l^d± ~ energy 
variation is: 



For the sake of simplicity, we limit our analysis to the isothermal case, = T. The 
entropy variation of the sub-layer is deduced from the two last equations, after introducing 
the conservation laws and after defining the polarized currents SI^ — {Isi — -^(ij.)/2, SIi — 
i^si - Idi)/'^, and the currents 7; = (7^ + /di)/2, I, = {I,^ + Isi)/2, 

^aa ^ pk-l^k _ pk^k+1 _ _^^k U/fc-l-fc _ + - 2*/) 

ditj 2 V / 

- \lA {is-'-"' - Is^'"-' - - I^h' {Sll-'-' - 611^'+' - 2^^, - 

- ^iil (/f-'-"*^ - + (30) 

where we have introduce the chemical potentials /ig = //^^ -|- //^j^, //J = /ig^ + /i^^, and 
the chemical affinities of the reactions, defined by A/Zg'^ = //^^ — //^j^ = ~'§^-i ^l^l' = 

,,k _ ,,k _ dE^ 

The entropy being an extensive variable, the total entropy variation of the system is 
obtained by summation over the layers 1 to ^2 where the layer 1 is in contact to the left 
reservoir and the layer Vt is in contact to the right reservoir W . 
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The total entropy variation is: 



dt ^' ^" 

n , Q 



k=2 k=2 



A:=2 



E ^(/^r' - /^t) 4"'^' + E ^/^«' + E M (31) 

fc=2 fc=l fc=l 

where the two first terms in the right hand side of the equality stand for the heat and 
chemical transfer from the reservoirs to the system E. 
The entropy variation takes the form 

= E^^^' + (32) 

i 

where Fi are generalized forces and X* are the conjugated generalized fluxes. The varia- 
tion of entropy is composed by an external entropy variation P^^*(t)/T and by an internal 
entropy variation dS^^^/dt. 

By applying the second law of thermodynamics dS^^^/dt > we are introducing the 
kinetic coefficients La/3 such that dS^"'^/dt = Yli^i (YIj ^ij^'^^ ■ By identification with the 
expression (j3ip . the kinetic equations are obtained, after performing the continuous limit, 



Jqs 




Jo I 








5Jf 




Ips 









Lss 

Lis 























-8 /Is 
dz 


Lu 
















—dAfis 
dz 





LssSs 












—d/iS 
dz 





Td 


Td 










dz 











T ^ 







A/is 














ri 







(33) 



The kinetic coefficients are state functions; L^- = Lij{S'' , N^, N'^) and the symmetrized 
matrix is positive : ^ {Lji + > 0. The coefficients Lint refer to the internal relax- 

ation processes j^J. According to Onsager relations, the kinetic coefficients are symmetric 
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or antisymmetric Lij = ±Lji. The coefficients are known from the two-channel model for the 
conductivity. The two last equations concern the internal (Lint) "density" variables ips and 
ipsd defined by = J^j. ip{z)dz. Due to the Curie principle, there is no coupling between 
spin polarized transport processes and the electronic transitions (the scalar process is not 
coupled to vectorial processes). 

VII. APPENDIX B 

This Appendix is devoted to the general resolution of the coupled diffusion equations 
Eqs. © in the case of a junction between two semi-infinite layers with the conditions of 
continuity written below: 



\ dz 



= 



(34) 



A/.f (0-) = A/xf (0+) 

Inserting the solutions given by © in Eqs. (jH^ and using Eqs. (0) the system becomes : 



6A+ + rfA_ = 



'^t(l+/3i) 



(lAL - 1/Al) + dXl - 
The h coefficient is given by : 



(35) 



h A_ 



A+ A_ 



1 1 



+ 



cr 



+ 



Jt 



(1 + /?,) ^ 



1 

A2 



iN 
''sf 



1 1 

a7^ aT 



(1 - m^y^iif^- 



(36) 



the d coefficient : 



d = Jt 



a^{l + P^)atA^ A_ 



(37) 
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and the a coefficient 



(38) 

For Isf ^ Isd '■ 

hfeJt jl + PsM-P^) .... 

1 /3.H(A^)^ ^^^^ 



For Z^/ > 



6^^Ul-/3f) ^ \.J,.,. (40) 
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